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Abstract 

We describe a non-perturbative method for computing the energy band structures of one- 
dimensional models with general point potentials sitting at equally spaced sites. This is done 
thanks to a Bethe ansatz approach and the method is applicable even when periodicity is 
broken, that is when Bloch's theorem is not valid any more. We derive the general equation 
governing the energy spectrum and illustrate its use in various situations. In particular, 
we get exact results for boundary effects. We also study non-perturbatively the effects 
of impurities in such systems. Finally, we discuss the possibility of including interactions 
between the particles of these systems. 
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Introduction 



The study of the so-called point or contact interactions covers an impressively large number of 
areas ranging from pure mathematics, where the main goal consists in defining and understanding 
them rigorously, to applied physics, where either they serve as good approximations for physical 
situations or they are created to control the behaviour of certain quantum systems. One particular 
case of such contact interactions but certainly the most well-known is the 5 potential which played 
a crucial role in two specific areas of special interest for us in this paper: solid state physics and 
exactly and integrable models (in one spatial dimension). Indeed, in the first context, it was used 
by Kronig and Penney £Q to model a free electron moving in a crystal lattice and is now a standard 
of textbooks in solid state physics. As for the second context, it has become one of the cornerstones 
of exactly solvable and integrable systems through the works of McGuire, Lieb and Liniger and 
Yang [21 El H] concerning a system of identical particles interacting through the 6 potential. In 
both cases, the corresponding quantum mechanical problem can actually be completely solved 
by taking a particular form of the eigenfunctions and one obtains the information on the energy 
spectrum. 

In the Kronig-Penney model, the periodicity of the potential and of the boundary conditions 
on the wavefunction allows to use the famous Bloch's theorem j^j which actually fixes the form 
of the eigenfunctions. In the Lieb-Liniger model, one assumes that the wavefunctions can be 
expanded on plane waves with appropriate coefficients that must be found. These approaches 
combined with the various conditions imposed on the wavefunctions of the problem lead to the 
allowed energy states of the models. Depending on the context, this gives rise to the famous 
energy band structures or to the famous Bethe ansatz equations [6 . 

Several generalizations of the two previous problems have been considered over the past 
decades. For instance, in Ej, the Kronig-Penney model is extended to the case where the 
6 potential is replaced by a more general point potential. In [Hj, boundary Bethe ansatz equations 
were derived by putting bosons with 6 interactions in a box while in [10] . impurity Bethe ansatz 
equations appeared by including a general external contact potential in a system of particles with 
5 interactions. 

With the advent of nanostructures and the ever increasing need for controlling quantum de- 
vices, the standard assumption of periodicity is no longer accurate enough. But then, Bloch's 
theorem no longer applies and, to the best of our knowledge, only approximate methods are used. 
In this paper, we mainly address the above problem and propose a general method to investigate 
the energy spectrum in one-dimensional models with equally spaced but otherwise arbitrary point 
potentials. Let us note that the breaking of periodicity can have various origins. For example, one 
can simply imagine that the periodic boundary conditions are replaced by Dirichlet or Neumann 
boundary conditions. Or, one could replace the S potential at one or more sites by a different 
point potential. We will see that the Bethe ansatz approach brings an elegant and more general 
alternative to Bloch's theorem in order to treat such problems. In all cases, we reduce the problem 
to finding roots of polynomials and solving Bethe ansatz equations. We also discuss the possibility 
of including interactions between bosons moving in such general potentials. 
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The paper is organised as follows. In section [TJ we discuss the general setup for treating free 
particles in an arbitrary equally spaced, external point potential. Then, in section EJ we validate 
our method in the context of periodic potentials by deriving a general result which encompasses 
previous results in the literature and in particular the well-known features of the Kronig-Penney 
model. Section |3] illustrates the use of our method in a context where periodicity is broken by 
the boundary conditions. This provides non-perturbative results for the boundary effects. This 
is further illustrated in Section 0] where we introduce impurities at one or several sites of the 
potential. Again the results are non-perturbative and this allows for a study of the effects of 
impurities with arbitrary strength. Finally, in Section we discuss the extension of our method 
to the case of bosons with 5 interactions moving in a general point potential. 



1 The general problem 

We study a one- dimensional system of free particles moving in the interval [—ML, ML], M ^ 1 
being an integer and L a length scale, with an external point potential sitting at each site x° = 
(M — 2j + 1)L, j = 1, . . . , M. It is sufficient to consider the one-particle Hamiltonian which takes 
the following form 1 



where, for j = 1, . . . , M 



Vj{X) 



Cj 8{x) + 4Aj ^0*0^ + ( 2 7i + Mfe) -^S(x) - (2 7i - i^) 5(x)^-, (1.2 



and Cj, Xj, jj, rjj 6 R. These very singular point potentials find their origin in the self-adjoint 



extensions of the free Hamiltonian — on the space Cq°(R\ {0}) of C°° functions with compact 
support separated from the origin. Written in this way, they are in fact quite formal but it is 
known that they are equivalent to imposing appropriate boundary conditions on the wavefunction 
at each site parametrized by U(2) matrices (see for example ^1] and references therein). 

In each region Rj : (M — 2j + ^4^)L < x < (M — 2j + ^-)L, the particle is actually free and 
we denote the wavefunction by <j)Ax) (see Figure |TJ). Now the idea is to impose all the boundary 
conditions (those corresponding to the potential at each x® and those at the ends of the interval) 
in a compact form by extending the approach of [10J - To do so, we collect all the pieces of the 
wavefunction into a single 2M-component vector defined for x e\ML — L, ML[ 

M 

$(ar) = <t>j{x + 2L{j - l))e J ®e+ + <j>j(-x+ 2L{M - j)) ej ® e_ , (1.3) 

3=1 

where {ej\l ^ j ^ M} is the canonical basis of C M and {e + , e_} is that of C 2 . We are now ready 
this paper, we use units such that H = 1 = 2m. 
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Figure 1: The nontrivial potential lies at the points — (M — — (M — 3)L, . . . , (M — 1)L (solid 
lines). The dashed lines correspond to a trivial potential used for convenience. 



to formulate the problem. As we said, this is just the free problem in each region Rj. In terms of 
$ this reads, 



dx 2 



$(x) = E®(x) , ML - L < x < ML. 



;i-4) 



Then, following for example [T2j, the M conditions corresponding to the general point potential 
are written as follows 



;i.5) 



(U + - I) $(x) + i(U + + I) ^-®(x) = for x^ML-L, 

ax 

where I is the 2M x 2M identity matrix and 

u + = 

\ U M J 

Each submatrix Uj ofU + corresponds to the potential Vj in (jl.lj) . It is a U(2) matrix parametrized 
by 



;i.6) 



£j G [0,7r), G C such that \fj,j\ 



1 . 



;i.7) 



where the symbol * stands for complex conjugation. The parameters in the matrix Uj are related 
to those of the potential (|1.2|) . For a discussion of the various parametrizations of point potentials 
see e.g. |Hj and references therein. Finally, the boundary conditions at the ends of the interval 
can be written in the same form as before 



(u- _ i) + i(u~ + I) —$(x) = for x — > ML 

ax 



where 



1 

1 



1 

1 



-e^ q is* 



1.9) 



^ ) 



The four coefficients at the corners of U also form a U(2) matrix 



u ' = s '{-k7;)' (L10) 

with the same constraints as in (|1.7p . They represent very general boundary conditions encom- 
passing the usual periodic, anti-periodic and box boundary conditions. These coefficients encode 
the behaviour of the wavefunction on the boundaries x = —ML and x = ML (see Fig. . The 

submatrices ( ^ J ^ of W _ correspond to the dashed lines in Figure [T] They are introduced for 

mere convenience in our approach. 

It is now obvious that one cannot use Bloch's theorem to solve this problem as explained in 
the introduction. The potential, albeit sitting on equally spaced sites, is certainly not periodic 
and we impose quite general boundary conditions. 

Instead, we formulate a Bethe ansatz 6 for $ 



$(x) = e ikx Ai + e ikx A R where 



M 

A i = E E A h e i * 

j=l e=± 
A/ 

^ = E E ^ j e * 



;i.n) 



In this form, is automatically solution of (jl.4j) with 

E = k 2 . (1.12) 

Inserting the ansatz in (jl.5j) and (jl.8j) . we find that $ is the eigenfunction we look for if and only 
if the 4M amplitudes A} ■ and satisfy the following relations 

^ = e 2ik{M - l)L Z + {-k)A n (1.13) 

^ = e 2ifcMi Z-(-A:)^ / , (1.14) 

where 

Z ± (k) = -[W ± - I - fc(Z/ ± + I)]" 1 [W ± - I + A;(W ± + I)] . (1.15) 

These 2M x 2M matrices characterize the type of potential and boundary conditions one is 
considering. Requiring a non trivial solution for the wavefunction, one finds the following equation 

det {Z + (-k) - Z-(-k)e 2ikL ) = . (1.16) 

This type of equation is usually called Bethe ansatz equations. To our knowledge, it is the first 
time that the Bethe equations has been established in this context. They impose constraints on the 
allowed values of k which incorporate the effect of the potential and of the boundary conditions. 
Solving in k as a function of the parameters controlling the potential and the boundary conditions 
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allows one to determine the energy spectrum and to study how one can modify it by tuning these 
parameters. As we will see, this has consequences on the energy band structure of the associated 
model. 

The above Bethe ansatz equations together with the method to get them for the very general 
model we are considering constitute the main result of this paper. They replace Bloch's theorem 
when periodicity is broken and provide a non-perturbative means to get the energy spectrum in 
such situations. 

2 Periodic potential 

Before exploring some applications of our method to more complicated cases in the following, 
we first show in this section that our approach consistently reproduce well-known results for 
periodic potentials and in particular for the Kronig-Penney model. We also obtain information 
on the bound states which are less studied. Finally, we compute the spectrum for an asymmetric 
potential which, to our knowledge, is not studied in the literature. 

2.1 General case 

In our language, a periodic potential with periodic boundary conditions is obtained by setting 



Ui 



and U„ 



■ = U 

1 

1 



M 



V 

-V [1* 



U 



(2.1) 
(2.2) 



In this case, the Bethe equations (jl.lfij) reduce to 



B+ 
T- 



T 
R 

— e 



— e 

2ikL 



2ikL 



T~ 



T+ 

R~ -e 2ikL 



— e 



2ikL 



Br 



—e 



2ikL 



where, for /jl = (j,r + ifii and v = ur + iuj, (fiR, fii, Ur, ujE 



I? 



Br 



(cos£ + fj,R)k 2 — 2i^ik — cos£ + fiR 
(cos ^ + fiR)k 2 + 2ik sin £ + cos £ — /xr ' 

(cos£ + fiR)k 2 + 2inik — cos£ + [j,r 
(cos£ + [iR)k 2 + 2ifcsin£ + cos£ — ur ' 



— e 
R 
T 



2ikL 



— e 



2ikL 



Br 



(2.3) 



-2vk 



(cos£ + [iRjk 2 + 2i/csin£ + cos£ — fiR ' 
2u*k 

(cos £ + ur) k 2 + 2i k sin £ + cos £ — jiR 



2 In comparison with the notation of [TO , R+ = R + (-k), T+ = T+(-k), R~ = R~(k) and T~ = T~(k) 



(2.4) 
(2.5) 
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The determinant can be seen as polynomial in e 2lkL of order 2M whose roots are functions of k 
of the form 

1 



Xz 



2 

2iir 



UJ 



*p 



T+ + lu p T~ ± ^J(u*p T+ - up T-) 2 + AR+R- 



p = 0,..,M- 1 



(2.6) 



where u = e "m is the M th root of unity. For M=l, we recover the result given in [TQ| . To get the 
spectrum, one has to solve in k the following equations 



„2ikL 



XI 



p = 0,..,M -1. 



(2.7) 



Let us introduce the shift operator Pf(x) = f(x + 2L). In the case under consideration, this 
operator commutes with the Hamiltonian and can therefore be diagonalized in the same basis. 
When acting on $, this operator is represented by the following 2M x 2M matrix 

/ 1 \ 



P 



\ 



(2.8) 



/ 



It is interesting to note that the eigenfunctions (jl.llj) constrained by relation ()2.7|) for a given p 
are actually common eigenfunctions of H and P as one can see from 



P $(x) = u p $(x) 



(2.9) 



This relation allows us to give a simple physical interpretation for p. Indeed, p is the usual crystal 
momentum 3 carried by the particle. In the usual approach, it labels the Bloch's functions and 
we recover here by our different approach that it is indeed a good quantum number for a very 
general periodic potential. We also remark that equation (|2.7jl is M-periodic in terms of p. This 
statement is equivalent to the standard restriction of the range of the crystal momentum to the 
first Brillouin zone. 



2.2 Kronig-Penney Model 

We now turn to the Kronig-Penney model to show that the equations (|2.7j) consistently reproduce 
the standard equations of this model. In this case, all the Vj's are given by a S potential with the 
same coupling constant c. This is obtained from (j2.1j) by taking 



U = -e* 
Therefore, we get 

T+ =T 



cos(f) isin(f) 
isin(f) cos(f) 



ik tan(£) 
1 + ifctan(f) 



^ G [0, 7r) and 



tan(£) = - . 

c 



and BP 



R- 



1 + ifctan(f) 



(2.10) 



(2.11) 



3 We choose the units so that H = 1 
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After some algebra, we see that the equations ()2.7|) become, for p=0,..,M-l and real k 



cos(2kL) 



fctan(£) 



sin(2A;L) = cos 



2np 



(2.12) 



These are the equations obtained by using Bloch's theorem (see for example [13]). They give the 
well-known plots that we show in Figure [2F- The plots represent the energy spectrum in terms 
of the crystal momentum p (restricted to the first Brillouin zone) for various values of £ (here 
M = 32 and L = 1). For £ = (circle on Figure EJ), the coupling constant, c, tends to +oo and 
the energy does not depend on p. The states are localized in the regions RjLi U Rj, j = 2, . . . , M 
which are separated by purely reflecting walls (there is no transmission). The choice £ = 7r/2 (box 
on Figure |2J) correspond to another particular case where the coupling constant vanishes. The 
states are delocalized and this is just the free particle model. The energy varies like the square 
of p. For an intermediate case, here £ = 0.2 (cross on Figure E]), the well-known energy band 
structure appears. 
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Figure 2: Lowest positive energy levels in 
terms of p for £ = O(circle); 0.2(cross); 
7r/2(box) 
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Figure 3: Bound states and lowest positive 
energy levels in terms of p for £ = 7r/2(box); 
2.3(cross); 2.5(diamond); 2.7(circle). 



Let us stress that in order to find equation (j2.12j) . we assumed k real. However, to study the 
possible bound states we have to look for purely imaginary solutions of ()2.7|) of the type k = iki, 
ki £ 1R. In this case, we are led to solve the following equations, for p = 0, .., M — 1 



1 + k I tan£)e 2klL = Man £ cos 



2np 



± Wl - fcf tan 2 £sin 2 



2ixp 



(2.13) 



Figure [3] shows the solutions of this equation (for £=2.3(cross); 2.5(diamond); 2.7(circle)) which 
gives negative energy levels. We represent also the lowest positive energy levels solution of (J2.12)) . 
The bound states appear only when £ > n/2 (i.e. c < 0). For £ — ► n~ , the coupling constant 
tends to — oo and the states are localized around the points x^. 



1 All the numerical resolutions and plots are realized with Maple 
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2.3 Asymmetric potential 

In this section, we study the periodic potential characterized by 



U 



e ia 
-e eia 



where e = — 1 (resp. e = 1) 



(2.14) 



This particular form of U for e = —1 (resp. e = +1) is obtained from (|2.1jl by setting // = 0, £ = 
and i/ = e ia (resp. /i = 0, i/ = 1 and £ = a). In this case, one gets 



(resp. i? H 



fc 2 -l 



k 2 + l 



R and T 



=F2A;e 



■iiia 



k 2 -l 



k 2 + 2ikt&na + 1 
For fc e M., equation ()2.7|) reduces to 



k 2 + 1 
iT and T ± 



^2A; 



cosa(/c 2 + 2iktana + 1) 



(2.15) 
(2.16) 



(A; 2 + l)sin(2fcL) = 2A;sin 



27TJ9 



a 



(resp. (k 2 + l)coso; sin(2A;L) + 2k sin a cos(2kL) = 2/csin 



2irp 
~M 



(2.17) 
(2.18) 



Figure |U (resp. Figure EJ) shows the behavior of the solutions of equation ()2.17j) (resp. equation 
()2.18j0 . Let us remark that, for a generic a, these potentials break the symmetry of the spectrum 
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Figure 4: Lowest positive energy levels in 
terms of p (e = —1) for a = (circle); 
7r/2 (cross); 7r(box) 



Figure 5: Lowest positive energy levels in 
terms of p (e = +1) for a = (circle); 
7r/4 (cross); 7r/2(box). 



under the transformation p — ► —p. We also remark that the behaviors of the two spectra are 
completely different. For e = —1, the gaps between the energy bands are independent of the 
parameters a whereas, for e = +1, they decrease when a increases. 
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3 Dirichlet or Neumann conditions 



In this section, we want to study the behavior of the energy band structure when we modify the 
boundary conditions. To illustrate this on a simple example we take the same potential as in (|2.1|) 
and simply change U q in ()2.2)1 to 



r I 1 



where e — — 1 



(resp. 



This implements the Dirichlet (resp. Neumann) boundary condition: 

0(0) = = (f)(L) (resp. 0'(O) = = <f>\L)) 



(3.1) 



(3.2) 



and breaks the periodicity of the model. This constitutes a basic example where our method still 
holds while Bloch's theorem breaks down. In this case the Bethe equations (jl.l6j) are equivalent 
to 



R + -ee 
T 



2ikL 



T 4 
R~ 



^2ikL 



-e 
R 
T 



2ikL 



T+ 

R~ -e 2ikL 



R- 



^2ikL 



-e 
R 
T 



2ikL 



R- 



ee 



2ikL 



0. 



(3.3) 



Again, this determinant is a polynomial of order 2M in e 2lkL whose roots read 



Y ± 



Y± = v / T+T r cos (^) ± JR+R--T+T- sin 2 (^) , 



-e(R + + R~) ± y/(R+ -R-) 2 + 4T+T- 



The Bethe equations governing the spectrum then read 

„2ikL _ Y ± ^ l^ q ^M-l 



2ikL 



q 

Y ± . 



(3.4) 
(3.5) 



(3.6) 
(3.7) 



Let us make a few remarks before showing the influence of the boundary on the energy spectrum 
for a particular example. First, the parameter q cannot be interpreted as the crystal momentum 
any more (the shift operator P does not commute with the Hamiltonian) . It does not label all 
the energy states since it is absent from equation ()3.7j) . On the other hand, only equation ()3.7|) 
depends on the parameter e characterizing the type of boundary conditions under consideration. 
We also stress that, although they look different, the roots in ()2.6|) and those in ()3.4|) can be 
written in the same form provided one replaces M by M/2 in ()3.4j) and relabels q appropriately. 
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From all this, we conclude that equations ()3.7|) give the energy states arising from the presence 
of the boundary while the rest of the spectrum will be identical to that of the periodic case. In 
particular, the boundary effects become negligible in the thermodynamic limit. 
Example 

We take once again the matrix U given by (|2.10j) to define the potential at each site but now, we 
impose Neumann boundary conditions. The Bethe equations ()3.6|) and ()3.7|) become respectively 
(for k real) 

cos(2kL) + — — sin(2A;L) = cos (—) , for 1 ^ q ^ M - 1 (3.8) 
/ctan4 \M / 

cos(2fcL) - A;tan(£) sin(2£;L) = 1 . (3.9) 

On Figure El below, we plot solutions of equations (|3.8|) and (|3.9|) for M = 16 and £ = 0.2. The 
solutions of ()3.8|) are represented by circles and those of ()3.9|) by boxes which we displayed at 
q = and q — 16 for convenience. To make the comparison even easier, we have also represented 
the energy band structure of the periodic case for M = 32 (crosses on Figure EJ) • We see that for 
the first and third energy band, the allowed energy levels coincide with the periodic case as the 
solutions of ()3.9j) nicely complete those of ()3.8j) . However, the effect of the boundary conditions 
shows up in the second energy band where two states are "missing" with respect to the periodic 
case. Again we recover graphically that this effect will become more and more negligible as M 
becomes larger. 
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Figure 6: Lowest positive energy levels in terms of q for £ = 0.2 
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4 Impurities 



So far, we have considered homogeneous potentials in the sense that the potential was the same at 
each site. As illustrated above, even with non-periodic boundary conditions, this does not deviate 
drastically from the periodic only local boundary effects arise. 

In this section, we want to explore very different situations which cannot be approached by 
perturbing the periodic case and therefore, show the advantages of our approach. Such situations 
occur when one modifies the potential at only one or several sites. Thus, these sites appear as 
what we call impurities. These models can be useful for describing localized defects in materials 
with one dimensional behaviour (nanowires, etc.) 

We illustrate this in the following by considering one such impurity and then a subset of equally 
spaced identical impurities in a given periodic potential. 



4.1 One impurity 



We restrict ourselves to the case of periodic boundary conditions and we imagine that the impurity 
is sitting at site 1 with a potential v± given as in (jl.2j) while all the other sites have the same 
potential v 2 — ■ ■ ■ — % = v. This model is obtained in our language by taking 



Mi 



Mi v\ 



U-2 



u 



M 



u 



fj, v 

-V* fX* 



and U„ 



1 

1 



(4.1) 



where U\ and U characterize respectively V\ and v. In this case, the Bethe equations (jl.lfij) reduce 
to 



Rt 



—e 



2ikL 



n 2ikL 



-e 
R 
T 



2ikL 



T 4 
R~ 



— e 



2ikL 



Rr 



2ikL 



n 2ikL 



n 2ikL 



-e 
R + 



Br 



(4.2) 



where Bf, T±, T{ and RT are given by relations (|2.4jl and (|2.5jl substituting ^, /i and v by ^i, n\ 
and ui, respectively. 

To simplify the expression of these Bethe equations, we consider the case of symmetric poten- 
tials i.e. R + = R~ = R, T + = T~ = T, Rf = BT = B± and T x + = T{~ — 7\. In this particular 
case, the vanishing of determinant ()4.2)1 is equivalent to the following equations 



2e 



2ikL 



y + 



T± 



T 2 + 4B 2 



(4.3) 
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where y is any root of one of the following polynomials 

T 2 (T X t Ri)(y 2M + 1) + T(R 2 -T 2 -R\- T 2 ± 2T 1 R)y M 

AI-l 

+R(2T 1 {R - Ri) ±(R- R^ 2 ± T 2 t T 2 ) ^ V 2j • ( 4 - 4 ) 

j=i 

Note that we consistently recover the periodic case when Ri = R and T\ = T since the polynomi- 
als ()4.4j) simply reduce to (y M — l) 2 . 

Example 

We imagine that both the impurity and the bulk potential are given by a 5 potential but with 
different coupling constants (i.e. the matrices U and U\ are given by (j2.10j) with parameter £ 
and £1, respectively). We study the behavior of the spectrum when the coupling constant of the 
impurity is varied while £ is kept fixed (we choose here £ = 0.2). 

Figures [7| and |H1 show the behavior of the energy levels of the first and second positive energy 
bands as functions of £i for M — 8. In the case where £i = £ = 0.2, we remark a degeneracy as 
the impurity becomes identical to the bulk potential. The corresponding energy levels coincide 
with those represented by crosses on Figure |2] Generally speaking, the effect of the impurity is far 
from being trivial. We see that in each band, two energy levels are severely modified depending 
on the strength of the impurity. This can have strong consequences on the corresponding system 
(if we think of conducting or insulating materials for instance). 




Figure 7: Energy levels of the first band in Figure 8: Energy levels of the second band 
terms of £i for £ = 0.2 in terms of £i for £ = 0.2 



4.2 TV identical impurities 

Here, we introduce periodically iV identical impurities (N G N) in a given periodic potential. From 
now on, let M be a multiple of iV i.e. M = MiN with Mi > 1. The impurities have potential V\ 
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U x = U Ml+1 = ■■■ = U (N _i )Ml+1 = ^ V \ ) (4.5) 



and the bulk potential is v. In this case, the matrices Uj take the following particular form 

U 2 = --- = U Ml = U Ml+2 = ■■■ = [/ (iV -i)M 1+ 2 = --- = U M = e i? ( ^ ) (4.6) 

As before, we restrict ourselves to the case of symmetric potentials for simplicity (R + = R~ = R, 
T + = T~ = T, Rf = R± = Ri and = = T x ). The Bethe equations take the form 



2e 2lkL = T (y + ± J T 2 - + 4R 2 . (4.7) 

Again, y is any root of one of the following polynomials, for ^ q ^ N — 1 

4^ P q + (y)P q -(y) + y 2Ml (7 9 - 1) V + l) 2 T 2 {R 2 -T 2 + T 2 - R\f , (4.8) 
where we have introduced 

lf{y) = T 2 (Ti =f i?!)(y 2Ml + 1) + i (V + 1) T(i? 2 -T 2 -R\- T 2 ± 21^)^ 

Afi-l 

+R(2T 1 {R - Rj) ±(R- R±) 2 ± T 2 t T 2 ) ^ y 2j , (4.9) 

i=i 

and 7 = e 2 ^" is the N th root of unity. The presence of 7 accounts for the inclusion of N identical 
impurities. In the case N — 1, one has 7 = 1 and the set of coupled equations ([4.7)1 . ([4. 8j) reduces 
to Q,g3J. 



5 AT bosons with contact interaction 

In this section, we extend the general method of Section [I] to M bosons with 5-interaction in 
presence of the potential defined in (jl.l)) . (jl.2)l . We restrict ourselves to the case of bosons for the 
sake of clarity but a treatment analogous to that in |H EI] > where the statistics is not given a 
priori, is possible. Our aim is to derive the corresponding Bethe ansatz equations which can be 
used then as above to get information on the energy spectrum of the model. However, we will 
not go into the details of exploring particular cases as the essential features of our approach have 
already been discussed in the previous sections. Instead, we trust that the interested people will 
adapt the general equations to their specific needs in order to study the influence of interactions 
in particular models with point potentials. 

We first need some definitions and notations. Let and 22J^- denote the permutation group 
and the Weyl group associated to the Lie algebra respectively. The group &jy consists of Af 
generators: the identity Id and M — 1 elements Ti, . . . , TV-i satisfying 

TjTj = Id , TjTi = T/Tj for \j - i\ > 1 , (5.1) 
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In particular, the last relation gives rise to the famous Yang-Baxter equation . For conve- 

nience, we denote a general transposition of by T^-, i < j, given by 

Tij = Tj-i . . . T i+1 TiT i+ i . . . Tj-i (5-3) 

The Weyl group 2Uv" contains 2- N 'M\ elements generated by Id, Ti, . . . , 2V-i an d -Ri satisfying 
(HHP . Q and 

Rx Ri = Id , (5.4) 

RyTxRyTi = T 1 R 1 T 1 R 1 , (5.5) 

= TjRy for j > 1 . (5.6) 

Let us define also Rj, j — 2, . . . , N as 

/,', 7} ; ...7; /,',•/; ...7} . . (5.7) 

Let xi,...,Xjf be the positions of the particles. The natural generalization of (|1.3|) leads to 
represent the wavefunction as follows 

4>(xx,...,xx) = cj)f 1 l[--'j^(xi,...,x Ar ) (5.8) 

in the regions, for 1 ^ k ^ M , 

Af-2j fc + ^-^L<x fc < ^M-2j fe + ^^^L, (5.9) 

and for xi,...,Xtf two by two different. Then, one defines the following (2M)^-component vector 
$(£!, ...,xjf)= Yl €u--Z (fjl^l • • • , f e h ® e £l ® ■ ■ • ® e iv ® e w (5.10) 

l<Jl,...,iAT<M 
ei,...,e^/-=± 

where, for 1 ^ k ^ A/", 

/^(x) = ±(x + 2kL-L- ML) -L + ML . (5.11) 

As before, the advantage of the vector <&(xi, . . . , xj^f) is that it contains all the information on the 
wavefunction for xy, . . . , x^ running only in the interval }ML — L, ML[ of length L. This allows 
to impose all the conditions for interactions between particles and with the external potential in 
a very compact form. 

Given a tensor product of spaces, (C 2 -^®-^, we define the action of a matrix A G End(C 2M ) 
on the k th space by 

A [k] = 1® •••<g>I®A<g>I<g> •••(8)1 . (5.12) 

fc-l N-k 
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We are now in position to write all the conditions the wavefunction of the problem has to 
satisfy: 

• Shrodinger equation: for Xj E\ML — L, ML[ and x\, . . . , Xjy two by two different 

M 

-^dl k <$>(x ll ... 1 x M ) = E<&(x u ... 1 x M ). (5.13) 

k=l 

• External point potential: for 1 ^ k ^ M 

[{U + - T) [k] + i(U + + l) [k] d Xk ] . . . , xjsf) = for ./,, — ML -L, (5.14) 

[(If- _ ][)[*] + + i)W aj . . . ; XN ) = o for i^ML. (5.15) 

• Interactions between particles: for Q e and 1 ^ i ^ A/" — 1, 

■ ■ .,x N )\ XQi = x+Q{w) =Q% Q~ l $(x 1; . . .,x N )\ XQi ^ Q(i+i) , (5.16) 
(d XQi - d XQ(i+1) ) . . . ,x N )\ XQi=x + (i+i) 

= Qf l Q- 1 [(d XQi - d XQ(i+1) ) + 2g] <&(*!, . . . , ar^)| XQi= *- (i+i) • (5.17) 

where Q is the usual representation of the element Q £ on (C 2M )®-^_ 

For bosons, the wavefunction <f>(xi, . . . , xjy) should be symmetric under the exchange of any two 
particles. In terms of $, this reads, for 1 ^ % < j ^ N, 

$(xi, . . . ,x u . . . ,Xj, . . .,x Af ) = T i:j $(xi, ...,Xj,...,x h .. .,xjj) ■ (5.18) 

Therefore, we can order the particles. In the following, we take ML — L < xi < ■ ■ ■ < x^f < ML. 
The ansatz consists in representing $ by 

$(zi,...,aw)= Yl e i{kpiXl+ - +kp " x " ) A P (5.19) 
Pe W M 



where 



A P = A pJuZm e n ® e ei <g> • • • ® e w <g> e eA . (5.20) 

ei,...,e^=± 

and, for any t> = (hi, . . . , vj^), 

v Ti = (vi,...,Vi-i,Vi +1 ,Vi,v i+ 2,...,vj^) (5.21) 

= (-vi,v 2 , ■ ■ ■ ■ (5.22) 

We need to determine the 2 J ^Af\(2M) J ^ parameters ^.p^'^v to find the solution of our prob- 
lem. Inserting in (jl.5|) and (jl.8|) . one gets 

A PRl = (Z+(-k P1 )) [l] e 2ik ^ M ^ L A P (5.23) 

A PRM = (Z-(-k PM )) W e 2ik ^ ML A P (5.24) 
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where Z ± is given by (jl.l5j) . From relations ()5.16|) . ()5.17|) . we get for 1 ^ j ^ M — 1 

Apt 3 = y{kpj - k P (j + x))Ap (5.25) 

where 

V(k) = ^. (5.26) 
k + ig 

Since our construction is based on 2U^, the Bethe ansatz is consistent if Z + and Y satisfy different 
relations (see for example ^UJ ) like the Yang-Baxter equation [H or the reflection equation 
[ISlEj. These relations hold true by direct computation with the explicit form of the matrices. 
The non-vanishing of the wavefunction implies the following constraints, for 1 ^ j ^ J\f, 



det 



. (5.27) 



(Z + (-%)) 111 - y{kj + h m )y{k 3 - k m )(Z-(-k,))^ l e 2 ^ 

These relations are the Bethe equations associated to our problem. We recover the results of 
[TUj by setting M — 1 and U q — ( ^ J J . As in the previous section, to get the spectrum of 
a given model with 5 interaction, one chooses the corresponding matrices Uj and U q and solves 

these equations in kj. The energy of the corresponding state is given by E = k 



2 
j' 

3=1 



Conclusions 

In this paper, we have presented a general method to address the question of the energy spectrum 
for a large class of one-dimensional models with equally spaced point potentials. We also describe 
how to extend the method to the case where particles interactions are present. The main results 
consist in Bethe ansatz equations which allow a non-perturbative treatment of these issues. 

We have illustrated the method in various typical situations of interest and have obtained exacts 
results on the effects of boundaries and impurities in such models. In view of the importance of 
the Kronig-Penney model in solid state physics, we believe that our method will serve as a useful 
toolbox to treat even more realistic situations which now occur in experimental situations with 
the advent of the " quantum technology" . 
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